arXiv: 1508.06797v2 [math.AP] 3 May 2016 


Lie Symmetry Analysis of the Black-Scholes-Merton Model for 
European Options with Stochastic Volatility 

A Paliathanasis* * 1 , K Krishnakumar 2 , KM Tamizhmanr 2 and PGL Leach^ 3,4,5 

1 Instituto de Ciencias Fi'sicas y Matematicas, Universidad Austral de Chile, Valdivia, Chile 
2 Department of Mathematics, Pondicherry University, Kalapet Puducherry 605 014, India 
3 Department of Mathematics and Institute of Systems Science, Durban University of 
Technology, PO Box 1334, Durban 4000, Republic of South Africa 
4 School of Mathematics, Statistics and Computer Science, University of KwaZulu-Natal, 
Private Bag X54001, Durban 4000, Republic of South Africa 
department of Mathematics and Statistics, University of Cyprus, Lefkosia 1678, Cyprus 

May 4, 2016 


Abstract 

We perform a classification of the Lie point symmetries for the Black-Scholes-Merton Model for European 
options with stochastic volatility, a, in which the last is defined by a stochastic differential equation with 
an Orstein Uhlenbeck term. In this model, the value of the option is given by a linear (1 + 2) evolution 
partial differential equation in which the price of the option depends upon two independent variables, the 
value of the underlying asset, S, and a new variable, y. We find that for arbitrary functional form of the 
volatility, a(y), the (1 + 2) evolution equation always admits two Lie point symmetries in addition to the 
automatic linear symmetry and the infinite number of solution symmetries. However, when a(y) = ao and 
as the price of the option depends upon the second Brownian motion in which the volatility is defined, the 
(1 + 2) evolution is not reduced to the Black-Scholes-Merton Equation, the model admits five Lie point 
symmetries in addition to the linear symmetry and the infinite number of solution symmetries. We apply 
the zeroth-order invariants of the Lie symmetries and we reduce the (1 + 2) evolution equation to a linear 
second-order ordinary differential equation. Finally, we study two models of special interest, the Heston 
model and the Stein-Stein model. 
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1 Introduction 

The Black-Scholes-Merton Model for European options is based upon some Ansatz for the stock price. Specif¬ 
ically, the process for the stock price is characterized by continuity, and it has the ability to hedge continuously 
with transaction costs and has constant volatility mmm- 
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In the Black-Scholes-Merton Model, the price of a financial asset is given by the soluton of the stochastic 
differential equation 

dSt = rStdt + aStdWt (1) 

where Wt is a Brownian motion, and the value u = u (t, S) of the option is given by the solution of the (1 + 1) 
evolution equation, 

-a 2 S 2 u y ss + rSu y s — ru + u y t = 0 (2) 

in which t is time, S is the current value of the underlying asset, for example a stock price, and r is the rate 
of return on a safe investment. The value of the option is subject to the satisfaction of the terminal condition, 
u (T, S) = U (S), when t = T. Finally, er is the volatility of the model. 

The Black-Scholes-Merton Model assumes constant volatility a. However, in real problems, er is not a 
constant. One possible generalisation of the model, Equation (0, is to consider that the volatility depends 
upon the time, t, and on the value of the stock, S, i.e., a = a (t, S). It has been proposed that a is a function 
of a mean Orstein Uhlenbeck process (3] . 

Consider that er = / (y), where y is given by the stochastic differential equation with the Orstein-Uhlenbeck 

term mm®-. 

dy t =a(m- y t ) dt + fidZ t (3) 

The new Brownian motion Z t can be correlated with Wt and be expressed as follows: 

Z t = pW t + \J 1 — p 2 Z t (4) 


in which Z t describes a Brownian motion independent of Wt and p is the correlation factor with values \p\ < 1. 

Hence, the Black-Scholes Equation © in the case of stochastic volatility is modified and the value u of the 
option is given by the (1 + 2 ) evolution equation 


(Mi + M 2 + M 3 + M 4 ) u (t, S,y) = 0 


(5) 


where the operators, Mi, M 2 , M 3 , M 4 , are dehned as follows: 


Mi 


I f 2 (v)S 2jf!_ + rS jL_ r+ <L 

2 ds 2 + b dS + dt 


M-2 


The function A (t, S, y) is 


rj2 O 

ppsf (y) , Ms = -/?A (t, S, y) —, and 

M> = \^b +a{m - v) h 


A (t, S, y) = p^rrr + 7 (t, S, y) \/l - p 2 

f{y) 


( 6 ) 

(7) 

( 8 ) 

(9) 


and u ( t , S, y) satisfies the terminal condition u ( T , S,y) = U ( S ) at time t = T. 

The operator Mi gives the Black-Scholes-Merton Equation © with volatility a = f (y ), M 2 expresses the 
correlation term between the two Brownian motions, Wt and Z t , of the European option and of the volatility, 
respectively, and M 4 is the Orstein Uhlenbeck process term. Finally, the term M 3 , the so called premium 
term, expresses the market price of the volatility risk [ 6 ]. The function 7 (t,S,y) in Equation © is the risk- 
premium factor which drives the volatility and follows from the second Brownian motion, Z t , where in the case 
of absolute correlation, i.e., \p\ = 1,7 it, S , y) does not play any role in the model. The first term of the rhs side 
of Equation © is called the excess return-to-risk ratio | 6 ]. The statistical importance of stochastic volatility 
has been confirmed in [ 8 .. 
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The purpose of this work is the study of the Black-Scholes-Merton Model with stochastic volatility, Equation 
©, by using the method of group invariant transformations, specifically the Lie (point) symmetries of the 
equation. The importance of Lie symmetries is that they provide a systematic method to facilitate the solution 
of differential equations because they provide first-order invariants which can be used to reduce the differential 
equations. Moreover, Lie symmetries can be used for the classification of differential equations. Furthermore, 
we can extract important information for the differential equation, consequently for the model, from the group 
of invariant transformations admitted. 

The first application of the Lie symmetries in financial modeling was performed by Gazizov & Ibragimov 
in ,9]. They studied the admitted group of invariant transformations for the Black-Scholes-Merton Equation 
0, with constant volatility and they proved that Equation © admits as Lie symmetries the elements of Lie 
algebra, {A 3j g © s A 3 p} © s 00 A 1 ( In the Mubarakzyanov Classification Scheme (121 IT, 12, T3J). This means 
that Equation © is maximally symmetric and according to the Theorem of Sophus Lie |14] there exists a 
transformation on the space of variables {t, S, u} in which Equation © can be written in the form of the 
heat equation. The last was an important result because the mathematical methods from physical science can 
be used for the study of differential equations in financial mathematics. A similar result has been found for 
the one-factor model of commodities HE which means that the three different equations, the heat equation, 
the Black-Scholes-Merton equation and the one-factor model of commodities equation, are equivalent at the 
mathematical level even if they describe different subjects. 

In recent years, Lie symmetries have covered a big range of applications in financial mathematics. For 
instance, the group invariants of the Cox-Ingersoll -Ross Pricing Equation have been studied in mi and the 
nonlinear Merton model in El- As far as concerns the Asian option, a Lie symmetry classification has been 
performed in [18]. As for generalisations of the Black-Scholes-Merton Model, the Lie symmetries and the 
reduction process of the nonautonomous model can be found in [19| l20j , while another generalisation of Equation 
© with a “source” was studied in [51] . 

Furthermore, in j55J [53] , the symmetry analysis of the space- and time-dependent one-factor model of 
commodities and of the nonautonomous two-dimensional Black-Scholes-Merton Equations were performed. 
For other applications of Lie symmetries in financial mathematics, see, for instance, [2411251 155] . and references 
therein. 

The stochastic volatility model, Equation ©, is a (1 + 2) evolution equation. Below, we perform a symmetry 
analysis and we determine the group invariant solutions. In particular, we restrict our analysis to the model 
in which the risk premium factor vanishes without necessarily \p\ = 1 and from Equation ©, only the term 
which expresses the return-to-risk ratio survives. Moreover, we study two models for European options with 
stochastic volatility, the Heston model m and the Stein-Stein model [28, • The latter is a model without 
correlation between the two Brownian motions, Wt and Z t , i.e., p = 0 in Equation ©. The plan of the paper 
is as follows. 

In Sectioi:© we give the basic properties and definitions for the Lie point symmetries of differential equations 
and we perform the symmetry classification for our model. We find that Equation © without the risk premium 
factor is always invariant under the {3Ai} 

oplus s ooAi Lie algebra. However, when / ( y ) is constant, Equation © is invariant under a larger Lie algebra. 
The application of the Lie symmetries to Equation © can be found in Section[3l in which we reduce the (1 + 2) 
evolution equation by using the zeroth-order invariants provided by the Lie symmetries and we derive invariant 
solutions. In Sections [4] and [5] we study two models of stochastic volatility for European options, the Heston 
model and the Stein-Stein model, respectively. For these two models, we find that both are invariant under the 
Lie algebra {3Ai} © s ooAi, and we apply the Lie symmetries to solve the equations of the two models. For the 
Heston model, the closed-form solution is expressed in terms of Kummer Functions, whereas for the Stein- Stein 
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model, the closed-form solution is expressed in terms of Hypergeometric Functions. Furthermore, we give some 
numerical solutions for the two models. Finally, in Section [Gl we discuss our results and draw our conclusions 


2 Lie Symmetry Analysis 


We consider the Black-Sclioles-Merton Equation with stochastic volatility governed by the evolution Equation 
([5j) for which the premium term depends only upon the return-to-risk ratio. For a time-independent rate-of- 
return, Equation ([5]) becomes 


H 


o = ^/ 2 (y) S 2 u,ss + PpSf {y) u S y + ^P 2 u,yy + 


+rSu t s + 


a(m- y) - /3p 


y — r 

f{y). 


u tV — ru + u t t 


Let $ be the map of an one-parameter point transformation such as 

( t,S,y )) = v! (■ t',S',y') 

with infinitesimal transformation (e is the parameter of smallness.) 

t' = t + (t, S, y, u) 


S' = S + e£ 2 (t,S,y,u ) 
y' = y + (t, S, y , u) 

u = u + ey (t, S, y, u) 


and generator 


dt' 

X = lT d < 

OS 


dS' 0 dy' 0 du' 0 

-d7 ds + ^ + ^ 


( 10 ) 

(11) 

( 12 ) 

(13) 

(14) 

(15) 

(16) 

(17) 


Consider now that u ( t , S , y) is a solution of Equation (fill) and under the map <F, Equation (fT2|) . v! (t', S' , y') 
is also a solution of Equation m ■ Then, we say that the generator X, of the infinitessimal transformation of 
the one-parameter point transformation, <f>, is a Lie (point) symmetry of ('Equation IIID and Equation (HID is 
invariant under the action of the map <f>. That means that there exists a function ■)/> such that the following 
condition holds [29] 

X® (U) = 4>H, mod (H) = 0 (18) 

or, equivalently, 

X [2] {%) = 0 (19) 

where X^ is the second prologation/extension of X in the space of variables {t, S , y , u, u t s,«, y , u y ss, u ,Sy> u yy}- 
Specifically X^ is defined from the following formula 


X [2] = X + rjf d u A + T)f-d u A 


where rjf, rj£j are given by the relations 


Vf = + rfr) A B ~ - u ?i u j£,B 

yfj = + 2 v A B(i u j) ~ £,% uA k + v A Bc u fi u j - 2 i%\B\ u f) uA 

- Z k BCU B iU C jU A k + V A B U% - 2 £%u$ )k - (u A k u% + 2ufjU^ )k J 


and 


( 20 ) 


( 21 ) 


( 22 ) 
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Table 1 : Lie brackets of the Lie 

symmetries of Equation (fill) for f(y ) 

= /o- 

X/,Xj] 

Xi 

x 2 

x 3 

x 4 

X5 

Xr 

0 

0 

—aX 3 

2/ 2 X 2 + (/ 2 - 2 r) X u 

olX 5 

x 2 

0 

0 

0 

2X U 

0 

X 3 

0X3 

0 

0 

0 

2af 0 X, 

X 4 

-2/ 0 2 X 2 - (/ 0 2 - 2r) X„ 

2X U 

0 

0 

0 

X 5 

—aX 5 

0 

-2af 0 X u 

0 

0 

X u 

0 

0 

0 

0 

0 


0 

0 

0 

0 

0 

0 


The importance of the existence of a Lie symmetry for a partial differential equation is that from the 
associated Lagrange’s system, 


dx l du A 
£* rj A 


(23) 


Zeroth-order invariants, ( x k ,u A ), can be determined which can be used to reduce the number of the 

independent variables of the differential equation. 

In the following, we perform a classification of the Lie symmetries of Equation m- Function / (y) is 
defined by the requirement that Equation m admit Lie symmetries. The latter requirement can be seen as a 
geometric selection rule as the Lie symmetries are generated from the elements of the Homothetic Algebra [30] 
of the (pseudo)Riemannian space, which defines the Laplace operator in the (1 + 2) evolution Equation (fill) . 
In our case, the (pseudo)Riemannian manifold is defined by the Brownian motions, Wj, Z tl of the stock price, 
S , and of the volatility, a , respectively. 

Before we proceed with the symmetry analysis, we remark that Equation E]) is a linear equation which 
means that it always admits the linear symmetry, X u = ud u and the inhnite-dinensional abelian subalgebra of 
solutions, Xt = b (f, S, y) d Ul where function b (t, S, y), is a solution of the original Equation (fill) l-IT . 


2.1 Classification 

From the symmetry condition Equation (1181) . we get a system of thirty-one equations ( for the derivation of 
the system, we used the symbolic package SYM of Mathematica 32, (33]) hr which the solution of the system 
gives the form of the generator Equation © of the transformation Equation m, that transforms solutions 
into solutions. From the latter system, we have the following results. 

For arbitrary function, / (y), Equation (1111) admits the Lie symmetries 

X 1= d t , X 2 = Sd s (24) 

plus the vector fields X u , Xf,. The algebra in which the Lie symmetries form is the{3Ai} © s ooAi. 

When / (y) = /o, Equation (fill) admits the Lie symmetries 

X-, =8 t , X 2 = Sd s , X 3 = e~ at d y (25) 

X 4 = 2f 0 tSds + 2^-fopdy + ((/ 0 2 - 2r) t + 2 In S) ud u (26) 

X 5 = e at (2 ftppSds + P 2 f 0 d y ) + (2af 0 (y - m) + 2/3 (p — r)) ud u (27) 

plus the vector fields X u , AV The Lie Brackets of the Lie algebra are given in Table [T] 
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We remark that the two-factor model of commodities is invariant under the same algebra of point trans¬ 
formations mm- That is an expected result because the two-factor model of commodities follows from the 
one-factor model in which the second factor, product, follows an Orstein-Uhlenbeck process. Moreover, as we 
discussed in the Introduction, the one-factor model is maximally symmetric just like the Black-Scholes-Merton 
Equation. 

On the other hand / ( y ) = f 0 means that the volatility er is constant. However, the second Brownian motion, 
Z t , in the space wherein a is defined, interacts with the Brownian motion Wt and modifies the Black-Scholes- 
Merton Model. However, in the case for which the correlation p vanishes, i.e., p = 0, Equation m is not 
reduced to Equation Q but only when the Orstein-Uhlenbeck process is identically zero, that is, /3 = 0, a = 0. 
Otherwise, the price u depends upon the Orstein-Uhlenbeck process. 

We continue with the reduction of Equation HID by applying the zeroth-order Lie invariants. Furthermore 
for every reduced equation we study the Lie symmetries. 

3 Group Invariant Solutions 

In this section, we apply the Lie symmetries in order to reduce Equation m- We study the two cases, 
f (y) = /o, and f (y) to be an arbitrary function. In order to perform the reduction and the later equation 
to give a solution of the original problem, there should be a constraint between the Lie symmetry vector and 
the terminal condition. However, we perform the reduction without considering the terminal condition at the 
moment because the initial conditions can be modified from different options. As far as the invariant solutions 
of the Black-Scholes-Merton Equation @ are concerned, see J3Tj. 

3.1 Arbitrary Function f (y) 

For an arbitrary functional form of / (y), as we saw above, Equation (1111) admits three Lie point symmetries in 
addition to the infinite number of solution symmetries. The last cannot be used for the reduction. Hence, we 
do not consider them. Moreover, a solution in which u does not depend upon one of the independent variables 
is not an acceptable solution, that is, the static solution is of no interest. Therefore, we perform reductions with 
the symmetry vectors Y\ = X\ + KiX u , Yi = X\ + KiX u and F 12 = X\ + cX 2 + n^Xu. 

Reduction with respect to the Lie invariants of the symmetry vector Y\ gives 

u(t,S,y)=exp[K 1 t]v(S,y) (28) 


where v (S, y) satisfies the equation 


0 = 7}f 2 {y)S 2 v,ss + Pl3Sf{y)v t sy + ^P 2 v,yy 


+rSv t s 


a{m- y) - /3p 


[x — r 

f (y). 


v >y - (r - «i) v 


(29) 


For this equation, we have that except the linear symmetry and the infinite number of solution symmetries 
(we call them trivial symmetries) the equation admits the vector field Y 2 = Sds, which is a reduced symmetry. 
Therefore, the application of Y 2 to Equation (l29l) gives the second-order ordinary differential equation 


P w ,vv 
where w = w (y) and 


28 

2a (to -y) + - (p« 2 / 2 (y) - p + r) 


, y + [(«2 ~ !) f 2 (y) + 2 ( rK 2 ~r + ki)] w = 0 (30) 

u {t, S, y) = S K2 exp [kR] w (y) (31) 
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Equation (1301) is a linear second-order differential equation, and it is well known that it is maximally sym¬ 
metric and is invariant under the special linear (si) algebra si (3 ,R) Lie algebra. 

Similarly, if we perform a reduction with Y 2 , the reduced equation admits the Lie Symmetries X\, X v , Xb , 
and finally the solution is again given by Equation m with the constraint Equation (l30l) 

Consider the application of the Lie symmetry vector Y 12 to Equation (1501) . We have that 

u (t, S, y) = exp [fc 3 t] v (z, y) , z = S'exp)—ct] (32) 


where 


0 


2 f(y) v, z 


+2 


ot(m — y) — pp 


2p/3f ( y) v, zy 

fj, — r 


f(y). 


+ P 2 V,yy +2 (r — 

v ,y - ( r ~ K 3 ) V 


c ) ZV t z 


(33) 


One can easily find that this equation only admits the Lie symmetry, zd Zl except the trivial symmetries, 
which is a reduced symmetry. Therefore, the application of the zeroth-order invariants of the symmetry vector 
( zd z + Kivd v ) in Equation (l33l) gives solution of the form Equation (IQTl) with the constraint Equation (l30l) 

We continue with the determination of the group invariant solutions for constant / ( y ). 


3.2 Constant Volatility 

For / ( y ) = /o, Equation (1111) admits six Lie point symmetries, plus the infinite number of solution symmetries. 
Moreover, Equation m is an (1 + 2 ) evolution equation, and, in order to reduce it to an ordinary differential 
equation, we have to apply the zeroth-order invariants of two Lie symmetries. From Table [L] we select reduc¬ 
ing Equation JEED by using the following two-dimensional subalgebras Aj = {Yi,Y 2 }, Bj = {Y 2 ,V 3 + kX u J, 
B h = {Y 2 , X 5 + kX u }, C t = {X 3 + kX u , X 4 } and C n = {X 4 , X 5 } . 

The reduction with the subalgebra Aj we studied in the previous subsection and the solution is Equation 
USD, where now from Equation (1301) we have 


w(y) 




(a (m-y) + f) 


+ w 2 A 


(_ C 2 _ 1 (a ( m-y) + f ) 2 
4a’2’ a 


We have that W\ and W 2 are constants, 

23 

Cl = — (pK 2 fo - p + r) , c 2 = [( k 2 - 1) fo + 2 (r-K 2 - r + ti 4 )] 
Jo 

and K, A are Kunnner Functions. 

We continue with the application of the remaining subalgebras. 

The application of Bi gives 

u (t , S , y) = S K2 exp ( ne at y ) ( t) 

where </>/ ( t ) is given by the first-order ordinary differential equation 


2fo4>i,t 


(/o «2 + 2 f 0 r) (k 2 - 1) + 2k/3 



p + r + 



JoP 2 k 


2 g 2 at 


4>i = 0 


(34) 


(35) 


with solution 


In 


<S>i (t) 

^/o 


(/o «2 + 2 r) 


(k 2 - 1 ) t - 


2k/3 

2 o/o 


( C 2 , atpf 0 \ at P 2 k 2 

r+ —) e r' 


(36) 
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From the subalgebra H//, we find the solution 


u {t , S, y) = pn {t ) S '* 2 exp ( —y 2 + 


(fte at - 2 am) 2 (r - p + pn 2 fl) 


P 2 h 


Pfo 


where pn {t ) is given by the expression 

ln ( ^//oO = [(/o « 2 + 2 ^)(k 2 - i) + 2 a] t 


+ 


ap 2 fl 


(P (r~p + pK 2 fo) + amfo) e~ 


2 at 


AafiP* 


For the subalgebra Cj , we have the invariant solution 

u (t, S, y ) = S+ (t) exp ( ne at y) pm {t ) 

where 




kPp 

af 0 t / 0 2 ' 2 ' 4/gt 


r 1 

72 + 2 


and 4>ni (t) is given by the expression 


ln 


( Pm = _ 
V Pmo 


1 


( 2p 2 — at) 


ln t + v r „ ' /3 2 K 2 e 2at 

a 2 t 


2 a/o 


Finally, from the subalgebra Cjj , we find the invariant solution 

Piv (t) 


u(t,S,y) = 


S^P’^ex p(V (t,y)) 


where 


V’ {t, V ) = 


^ (*, Z/) = 


p 2 {at-2p 2 ) V ° 

and function piv (t) is given by the expression 

Piv (t) 


sjat — 2 p 2 

4 pf 0 a {m-y)+P(a (/ 2 -2r)t- 4p {p - r)) + aP 2 
2pf 2 {at - 2p 2 ) 

2/3 (r - p + f -fp - 2 rp'j + 2/ 0 am 


P 2 fo {at-2p 2 ) 


ty 


ln 


Piv o 


fi (8a - 4 r - /J) - 4r> (<rf2 _ y() + (/o’ - 2-)' 


8/o (a* - 2p 2 ) 


2fo a {at — 2p 2 ) 


2 (/o- 2r2 ) {foam - P {p - r)) 3 (4/3 (r - /i) + 2/ 0 am) 2 

H- - ———:-- P ~ 


flap {at - 2p 2 ) 


P 2 fo {at - 2p 2 ) 


-mp 


2 (m - H 

fo a {at — 2p 2 ) 1 


(37) 


(38) 


(39) 


^2amf + 2p (r — p — 2pr + ^/o) + “—( 40 ) 


(41) 

(42) 

(43) 


(44) 


In the following section, we study a special model for stochastic volatility which has been proposed by Heston 


m- 




























4 Heston Model 


In the Heston model for stochastic volatility the stock price, S , and the volatility, a = f (y), satisfy the stochastic 
differential equation given below 


dSt = rSt.dt + St'jYtdWt and, 
dY t = 9 {fh — Y t ) dt + d\fY t dZ t 


(45) 

(46) 


where, in comparison with Equations © and m, we observe that / (y) = y/y and /3 = 6^/y. The differential 
equation which corresponds to that model is 


0 = ^YS 2 u,ss + pSYSu,sy + ^S 2 Yu : yy + rSu,s 

+ (9 (fh — Y) — XY) u t Y — ru + u tt 


(47) 


Before we proceed with the symmetry analysis of Equation (S3, we perform the coordinate transformation 
y = y 2 . Then, Equation (1-T71) becomes 


0 = 


\y 2 S 2 u tS s + PpySu,sy + ^/3 2 u, V y 


rSu 


,s 


1 

+ 2 


, i C2 

ay+ — 
y 


u tV — ru + u t t 


in which we have made the replacements 


/3, ci = (9 — A) and C 2 = 9m — 


(48) 

Equation C51) can be 

compared with Equation © for / ( y) = y. However, the risk premium factor of Equation © is not zero and 
has absorbed the term of the Orstein-Uhlenbeck process. 

From the Lie symmetry condition, Equations (11811 for (14811 we find that this equation admits the Lie symmetries 


X\ — dt , X2 — Sds 


(49) 


plus the X u , Xb, that is, Equation (l48l) is invariant under the Lie algebra {3Hi} © s 00 H 1 . 

We continue with the application of the Lie symmetries in order to reduce Equation (147)1) . We follow the 
results of Section [3© that is, we apply the group invariants of the subalgebra {Fi, Y^}- 
We find that the corresponding invariant solution of Equation (l48l) is 

u(t,S,y) = S K2 e^ t W(y) (50) 

where W (y) satisfies the linear second-order differential equation: 

P 2 W, VV + ^2ac 2 /3 p + ci y + jj W iV + (2/ci - 2 r (1 - k 2 ) + y 2 {k 2 2 - «?)) W = 0 (51) 

the solution in closed form of which is expressed in terms of Kummer Functions. 

In Figures [Hand© we give numerical solutions of Equation (1511) . Figure[l]is for negative value of Aft, whereas 
Figure© is for positive value of Aft, where Aft = k 2 — sq. 

5 Stein-Stein Model 

The model which has been proposed by Elias M. Stein and Jeremy C. Stein [28] describes an European option 
with stochastic volatility for which the correlation among the two Brownian motions vanishes, i.e., p = 0, in 
Equation ©. Moreover, they considered that the risk premium factor is constant, i.e., j(t,S,y) = 70 and the 
volatility is / (y) = y, while the stochastic differential equation is 

dSt = rStdt + StYtdWt and 
dY t = 9 (fh — Y t ) dt + SdZ t 
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(52) 

(53) 




Evolution of W(y) for p>0, <^<0 and c 2 <0 


Evolution of W(y) for p>0, c^O and c 2 >0 



Figure 1: Evolution of the solution Equation (l5lT) of the Heston model. For the numerical solutions, we select 
p = 0.5, /3 = |0.7|, r = 0.5, K\ = 1 and K 2 = 0.5. The left figures are for negative ci, C 2 , while the right figures 
are for negative c\ and positive C 2 - The solid lines are for |c 2 1 = 5 |ci| , the dotted lines are for |c 2 1 = 0.2 |ci| , 
and the dash-dash lines are for |C 2 1 = |ci|. The top figures are for /3 > 0, while the lower figures for /? < 0. 



Figure 2: Evolution of the solution Equation (ISlT) of the Heston model. For the numerical solutions, we select 
p = 0.5, /3 = |0.7|, r = 0.5, K\ = 1 and «2 = 1.01. The left figures are for negative ci, C 2 , while the right figures 
are for negative c\ and positive C 2 . The solid lines are for |c 2 1 = § |ci|, the dotted lines are for |c 2 1 = | |ci| , 
and the dash-dash lines are for |C 2 1 = |ci|. The top figures are for /? > 0, while the lower figures for /? < 0. 
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Figure 3: Evolution of the solution Equation (1571) of the Stein-Stein model. For the numerical solutions, we select 
/3 = 0.5, r = 0.5, K\ = 1 and k 2 = 0.5. The figures are for ui = 0, to = 0.5, u> = —0.5 and u> = a, respectively. 
The solid lines are for a = 0.1, the dotted lines are for a = 0.3 and the dash-dash lines are for a = 0.4. 


Therefore, from Equation d5j) , we have that the (1 + 2) evolution differential equation of the Stein-Stein 
model is 

-y 2 S 2 u t ss + -/3 2 u y!/ + rSu t s + [a{m-y) - /Sqo] v )V -ru + u >t = 0 (54) 

From the Lie symmetry condition Equation (1181) . we find that Equation (1541) admits the Lie symmetries 


X! = dt, X 2 = Sd s , X u = ud u , X b = b (t, S , y ) d u (55) 

which form the Lie algebra {3Ai} © s 00 .A 1 , and it is the admitted algebra of the Heston model and of Equation 

m for the arbitrary function / (y). 

Following the steps of the previous sections, we find that the invariant solution of the Stein Stein model 
with respect to the Lie algebra {Yi,!^} is 

u {t, S, y) = S K2 e Klt Y (y) (56) 

where Y (y) is given by the linear second-order differential equation 

P 2 Y )VV + (2 u - 2 ay) Y y + (2r (k 2 - 1) + 2m + y 2 - ref)) Y = 0 (57) 

The closed-form solution of this equation can be expressed in terms of the Hypergeometric Functions, where 

ui = am — (3 7 o- In Figures [3] and 01 we give the numerical evolution of Y (y) for various values of the parameters, 
w and a, for negative Ak and positive Ak, respectively, where A k = k 2 — Ki- 

6 Conclusions 


Volatility with a stochastic process has been shown to be essential for Financial Mathematics. In this work, 
we studied the algebraic properties, i.e., the Lie symmetries, of the modified Black-Scholes-Merton Equation 
for European options with a stochastic volatility. We have shown that the autonomous model without the risk 


11 

















Evolution of Y(y) forto=0 


Evolution of Y(y) forw>0 
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Figure 4: Evolution of the solution Equation (l57l) of the Stein-Stein model. For the numerical solutions, we 
select /3 = 0.5, r = 0.5, = 1 and K 2 = 1.01. The figures are for w = 0, u = 0.5, ui = —0.5 and w = a, 

respectively. The solid lines are for a = 1.1, the dotted lines are for a = 1.3 and the dash-dash lines are for 
a = 1.4. 

premium factor is invariant under a group of point transformations which form the {3^4i} ® s ooAi Lie Algebra 
for an arbitrary functional form of the volatility, a. Moreover, when the volatility is constant but the price 
of the option depends on the second Brownian motion, in which the volatility is defined, the modified Black- 
Scholes-Merton Model is invariant under six, plus the infinity, Lie point symmetries and it is not maximally 
symmetric as the Black-Scholes-Merton Equation with nonstochastic volatility is. 

Furthermore, we showed that the Black-Scholes-Merton Model, in which the volatility is constant, but is 
defined by an Orstein-Uhlenbeck process, is invariant under the same group of point transformations as that 
of the two-factor model of commodities. The reason for that is that the two models have in common the terms 
which follow from the Orstein- Uhlenbeck process. 

Moreover, we applied the zeroth-order invariants of the Lie symmetries, and we reduced the model to a 
linear second-order differential equation. As far as the case of constant volatility is concerned, we found the 
closed forms of the group invariant solutions. 

Finally, we studied the algebraic properties and the invariant solutions of two models, the Heston model 
and the Stein-Stein model, with stochastic volatility of special interest. For each model, we found the invariant 
solution and we gave some figures for the evolution of the models. Of course because Equation © is a linear 
equation, the general solution is given by the linear combination of the invariant solutions that we have found, 
while the latter are constrained by the initial conditions and the boundary conditions of the model. 

A general consideration of Equation (f5|) . in which the risk premium factor plays a role is still in progress, 
and the results will be published elsewhere. 
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